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ABSTRACT 

We compare predictions of large-scale cosmological hydrodynamical simulations 
for neutral hydrogen absorption signatures in the vicinity of Mq haloes 

with observational measurements. Two different hydrodynamical techniques and a 
variety of prescriptions for gas removal in high density regions are examined. Star for¬ 
mation and wind feedback play only secondary roles in the Hi absorption signatures 
outside the virial radius, but play important roles within. Accordingly, we identify 
three distinct gaseous regions around a halo: the virialized region, the mesogalactic 
medium outside the virial radius arising from the extended haloes of galaxies out to 
about two turnaround radii, and the intergalactic medium beyond. Predictions for the 
amount of absorption from the mesogalactic and intergalactic media are robust across 
different methodologies, and the predictions agree with the amount of absorption ob¬ 
served around star-forming galaxies and QSO host galaxies. Recovering the measured 
amount of absorption within the virialized region, however, requires either a higher 
dynamic range in the simulations, additional physics, or both. 

Key words: cosmology: large-scale structure of Universe - quasars: absorption lines 
- galaxies: formation - intergalactic medium 


1 INTRODUCTION 

The environment of a forming galaxy is expected to be vi¬ 
olent and complex. Cooling and accreting gas will be ther¬ 
mally unstable on multiple scales, leading to star formation. 
Supernovae-driven winds from young massive stars will drive 
back accreting gas. An active galactic nucleus (AGN), par¬ 
ticularly in more massive galaxies, will produce jets that will 
further drive out gas, while the accompanying compression 
of gas latitudinal to the jet may enhance star formation. The 
evolution will be occasionally punctuated by major mergers 
with other haloes. 

Cosmological simulations have begun filling in the 
physical details of these processes. A picture is emerg¬ 
ing of two modes of gas accretion, a “hot mode” in 
galaxies sufficiently massive for most gas to shock heat 
to virial temperatures before the onset of gradual cool¬ 
ing, and a “cold mode” in galaxies in which much 
of the gas does not sustain high post-shock tempera- 
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tures and instead flows in primarily as cool gas streams 
j[Kjy_ef^jjj|2000| jKei^^ akEoPSl : iDekel k, Birnbolmll2006l : 


Ivan de Voort et al.l 201l|), ascenario sup ported by theoreti¬ 
cal arguments i Birnboim fc Dekell 


The cooling gas feeds star formation in the galaxies. 
Simulations, however, tend to over-predict the amount of 
stars formed. Supernovae have long been suspected of sup¬ 
press i ng star formation in galaxies, from ellipticals (iLarsonl 
I 1974 I : I Silk et all 1 19861 ) to dwarfs llPekel fc SilkI Il986h . to 
bring the amount of stars formed into agreement with ob¬ 
servations. This view has been backed by large-scale hydro¬ 
dynamical simulations, strongly suggesting that supernovae- 
driven wind f eedback may regulate the star formation rate in 
galax ies (e.g. lKav et al.ll2002 : ISchave et ani2010l : lHaas et ^ 
I2OI3I) . AGN jets may furthermor e regulate star formation 
in the more mas sive galaxies (e.g. Ivan de Voort et al.ll201ll : 
iHaas et"^ lliH). 


However, direct observational evidence for these mech¬ 
anisms is scanty. While hydrogen and metal emission lines 
show evidence for winds in star-forming galaxies at red- 
shifts z > 1 with velocities excee ding the escape velocities 
of the haloes llGenzel et al.|[20ll] ). the extent of the winds 
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and their impact on the surrounding gas are still unknown. 
Maps of gas emission in active galaxies reveal winds extend- 
ing over tens of kilopa rsecs, but so far only at low redshifts 
dHarrison et al.l l2014l l. Consequently, the search for cold 
streams in absorption in the spectra of background Quasi- 
stellar Objects (QSOs) or bright galaxies has attracted in¬ 
creasing recent attention. The availa bility of large QSO cat¬ 
alogs with a hig h sky density (e.g. lAbazaiian et al.l l2009l : 
lAhn et ^ l2012|j . has presented an opportunity to study 


the gaseous environments of star-forming galaxies (SFGs) 
and AGNs through absorption measurements with unprece¬ 
dented detail. Observations of intergalactic Hi show a rise 
in the Hi optical d e pth near Lyman-break galaxies (LBGs ) 
dSteidel et ah 20]T|;^Crjght^ et_M.l l20lI] : lRakic et al.ll2012l i 
and QSOs ( Prochaska et al.ll2013l l. and thus evidence for 
extended cool gas surrounding galactic haloes. 

Several galaxy formation simulations have sought to 
establish consistency with these data. Of particular fo¬ 
cus have been absorbers optically thick at the photo¬ 
electric edge, Lyman limit systems (LLS) and damped 
Lya absorbers (DLAs). Simulations suggest a large frac¬ 
tion of these absorption systems at 2 > 2 may be 

acco unted for by cold-mode accre t ion of circumgalactic 
gas dFau cher-Giguere fc Kered 20ll]: iFu magalli et ^l201ll : 


Ivan de Voort et^ 20121 : Shen et 2013l L However, nu¬ 

merical models have been less successful at reproducing the 
observed absorption around QSOs relat ive to SFGs (but 
see recent work by iRahmati et ^ l2015l L While the cov¬ 
ering fraction of optically thick systems of SFGs in sim¬ 
ulations matches the observations within the errors, most 
previous simulation predictions fall short of the covering 
fractions measured around QSOs dFumagalli et al.l I2OI4I : 
iFaucher-Giguere et al.|[2014l L Agreement with the trend of 
increasing Hi absorption with decreasi ng impact parameter 
is also reported for the SFG sample of |Steidel_et_alJ (Hoiol) 


dFumagalli et al.ll201ll : iGoerdt et al.ll2012l : Shen et ^l2013ll , 

but for a sample of QSOs, simulations greatly underpredict 
the observed amount of absorption by circumgalactic gas 
dProchaska et al]l2013l 'l. 

In this paper, we examine the Hi absorption around 
galaxy haloes in large scale cosmological hydrodynamical 
simulations, by considering both the circumgalactic gas and 
the extended region surrounding it as it merges into the in¬ 
tergalactic medium (IGM). Previous simulations of the cir¬ 
cumgalactic medium (GGM), generally defined as gas within 
300 kpc (proper) of a galaxy, have instead typically been 
zoom-in models focussed on the immediate surroundings of a 
galaxy, extending out to at most a few virial radii. They have 
not therefore been able to fully exploit the growing amount 
of data probing the gas around galaxies to distances of 
severa l comoving megaparsecs. Exceptions are iRakic et al.l 
ll2013tl . who used the medi an Lya optical depth measure¬ 
ments of iRakic et al.l (l2012tl to constrain the masses of the 
SFG haloes in the K eck Baryonic St r uctur e Survey, and a 
very recent studvJjv jRjahmjti el^l.1 (l2015h using the EA¬ 
GLE simulation JSchave et al. 2015lj . 

The analysis presented here is inst ead based on the Enzo 
and GADGET-3 simulations described in iMeiksin et al.l (l2014l l 
(hereafter Paper I). These are of sufficiently high resolution 
and volume to recover the hydrogen absorption properties 
of the IGM. In Paper I we demonstrated that these IGM 
simulations make consistent predictions for the intergalactic 


gas properties beyond the halo turn-around radii, despite 
their differing numerical algorithms. However, the gas prop¬ 
erties on smaller scales are highly dependent o n star forma¬ 
tion and feedback impl ementations (see also ISuresh et al.l 
I2OI5I : [Nelson et '^l 2015 l L We shall use these simulations to 
construct spectra to examine how well the models match 
the observed statistics of Hi absorption around galaxies and 
QSOs. We differ from previous work, however, by first taking 
a step back and asking whether the absorption signatures of 
gas around galaxies show any evidence for star formation 
even without feedback. We then seek to establish the physi¬ 
cal extent of excess hydrogen absorption over the contribu¬ 
tion from the diffuse IGM around galaxies. We finally ask 
what observational evidence hydrogen absorption signatures 
provide for wind feedback in addition to star formation. 

This paper is organized as follows. In Sec. [2] we sum¬ 
marise the numerical simulations used in this work. The 
resulting Hi optical depths and covering fractions around 
haloes are presented in Sec. [31 and these results are dis¬ 
cussed in detail in Sec. (4) We present our conclusions in 
Sec. [SI and simulation convergence tests are discussed in an 
Appendix. All results are for a flat AGDM universe with 
the cosmological parameters flm = 0.28, Qbh^ = 0.0225 
and h = Ho/100 km s“^ Mpc“^ = 0.70, representing the 
present-day total mass density, baryon density and Hub¬ 
ble constant, respectively. The power spectrum has spectral 
index Us = 0.96, and is normalized to ag = 0.82, consis¬ 
tent with the 9- year Wilkinson Micro wave Anisotropy Probe 
(WMAP) data ( Hinshaw et al.ir2013h . 


2 NUMERICAL SIMULATIONS 

2.1 Cosmological hydrodynamics codes 

The numerical sim ulations are performed using ENZO, v.2 
d Bryan et aLllioi^ L and GADGET-3, an updated version of 
the publicly available code GADGET-2 (last described by 
ISpringel|l2005l 'l. The runs are summarized in Table [T] and 
further details may also be found in Paper 1. As we fo¬ 
cus on high redshift haloes in this work, all the runs were 
stopped at 2 = 2. The Eulerian ENZO simulations were per¬ 
formed on the top-level grid only. One simulation uses 512® 
mesh zones and 512® cold dark matter particles in a box 
30 Mpc (comoving) on a side (E30_512), and a second has 
1024® mesh zones and 1024® cold dark matter particles in 
a 60 Mpc box (E60_1024). The grid zones have a comov¬ 
ing width of 58.6 kpc. Since ENZO uses the particle-mesh 
method, the force resolution is limited to two grid zones, 
adequate for resolving the Jeans length of the photoionized 
gas. The smoothed particle hydrodynamics simulations per¬ 
formed with GADGET-3 were run with 512® cold dark matter 
particles and 512® gas particles in a box 30 Mpc (comov¬ 
ing) along a side, with a force softening scale of 1.4 kpc 
(comoving). The baryonic and dark matter particle masses 
are 1.3 x 10® Mq and 6.4 x 10® Mq, respectively. 

No gas is removed from the ENZO simulations, produc¬ 
ing a “no star formation” control sample. However, be¬ 
cause of the high densities and temperatures reached in 
the Lagrangian GADGET-3 simulations, gas must be removed 
in rapidly coo ling regions. We adopt the star formation 
prescription of ISoringel fc Hernauistf (I 2 OO 3 II (SH03) in one 
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Table 1. Summary of the simulations performed in this work. The columns, from left to right, list the simulation name, the box size in 
comoving Mpc, the number of resolution elements in the simulation, the code used for the run, the star formation prescription and whether 
or not the model includes supernovae-driven winds. 


Name 

Box size 

Resolution 

Method 

Star formation 

Winds 


[Mpc] 

elements 




GSOqLyo 

30 

2x512^ 

GADGET-3 

qLya 

N 

GSOsfnw 

30 

2x5123 

GADGET-3 

SH03 

N 

GSOsfw 

30 

2x5123 

GADGET-3 

SH03 

Y 

E30.512 

30 

5123 

ENZO 

none 

N 

E60.1024 

60 

10243 

ENZO 

none 

N 


model (GSOsfnw), and include their supernovae feedback 
prescription with a wind velocity of 484 km in another 
(GSOsfw). Note, however, that any conclusions drawn from 
its use are exclusive to this part icular algorithm (see also 
e.g. lDalla Vecchia fc Schavell200^ ~). Indeed, absorption mea¬ 
surements like those we examine may provide a basis for 
discriminating between different feedback models. A com¬ 
parison simulation, discussed in the Appendix, is also per¬ 
formed using the “quick Lya” method (GSOqLya), for which 
all gas cooler than 10® K and with an overdensity exceeding 
1000 is converted into collisionless “star” particles (without 
feedback). This is a computational trick used in IGM simu¬ 
lation s to significantly speed up the computation dViel et al.l 
l2004h . We use it as an alternative for testing the sensitivity 
of the results on the gas removal scheme. 

Both the ENZO and GADGET-3 codes solve the time- 
dependent ionization equations for hydrogen and helium. 
Photoionization, collisional ionization and radiative and 
dielectric recombinations are included. Thermal balance 
includes cooling by radiative and dielectronic recombi¬ 
nation, collisional ionization and excitation and thermal 
bremsstrahlung losses as well as Compton cooling off the 
cosmic microwave background. The cod es adopt the pho¬ 
toioni zation and photoheating rates of iHaardt fc Madaul 
(l 2012 ll . except we adjusted the photoionization heating rate 
of singly ionized heli um to reproduce th e evolution of the 
IGM temperature in iBecker et al.l ll201lll for 7 = 1.3. The 
codes used identical at omic rates for the heating and cool¬ 
ing, as summarized in iMeiksinI (l2009fl . except for adopting 
the Hi electron excitat ion and collisional cooling rate of 
IScholz fc Waited (Il99lh . We have not included any pho¬ 
toionization from the central galaxies or QSOs, which could 
reduce t he amount of absorption from high c olumn density 
systems (iMiralda-Escudd I2OO5I : ISchavel I2OO6II . However, to 
date there is no firm evidence for a widesprea d transverse 
proximity effect fr om either SFGs or QSOs (ICro ftl |2004 
iKirkman fc Tvtl^ 12008(1 . We also do not include the self¬ 
shielding of Hi absorbers in our main analysis, although in 
the Appendix we show this has only a small effect on the 
integrated absorption statistics examined here. 


2.2 Extraction of Lya absorption statistics 

The Lya absorption spectra are computed at restframe ve¬ 
locity V along a line of sight through a periodic box of co¬ 
moving side L as exp[—ra(i>)], where 


Ta{v) = ^ f dxnHiix)b{x) . ( 1 ) 

Here uios = xH{z)L/{l + z) + Vpec at dimensionless dis¬ 
tance X, in box side units, for Hubble parameter H{z) 
at redshift z and line-of-sight peculiar velocity Upec The 
Doppler parameter b{x) = [2kBT{x)/m-iiY^‘^ for gas at 
temperature T{x), CTq, = {ne^lmec)fa is the total cross- 
section for Lya scattering with upward oscillator strength 
fa — 0.4162, Ac — 1215.67 A is the Lya wavelength, 
fcs is the Boltzmann constant and mn is the mass of 
a hydrogen atom. For the ENZO simulations, the gridded 
Hi density, temperature and peculiar velocity fields are 
used. For the GADGET-3 simulations, the spectra are com¬ 
puted along rays using H i-weight ed contributions from 
the particles (iHernauist et al.lll996l l. To match the resolu¬ 
tion of the observations, we smooth the resulting spectra 
with a Gau s sian o f FWHM 8 km s“^ for comparison with 
iRakic et al.l (l20 12ll an d FWHM 125 km s“^ for comparison 
with IProchaska et al.l (l2013l f . 

We extract spectra for a range of impact parameters 
around haloes identified in the simulations using the grid- 
based halo finding algorithm described in Paper I. We con¬ 
sider three spectral signatures and their dependence on im¬ 
pact parameter: the 2D absorption as measured by the me¬ 
dian Lya optical depth around the haloes, the Lya absorp¬ 
tion equivalent width and the fractional deviation of the Lya 
absorption from the mean intergalactic value. Note the lat¬ 
ter two statistics are integrated quantities. The first of these 
is the equivalent width around a halo at measured velocity 
Uhaio along a line of sight to a background QSO, which is 
computed as 





ui(fex. An) = — / dv \ \ 

(2) 


J 



over a velocity window of width Av centred on the position 
of the halo, displaced tra nsversely by an amount b±. Fol¬ 
lowing |Rakic|eF^ (l 2012 fl . we also consider an alternative 
determination of the equivalent width by dividing the trans¬ 
mission by a factor meant to correct for the errors incurred 
in continuum fitting low-resolution spectra. The second in¬ 
tegrated absorption statistic we consider is the fractional 
deviation, of the a bsorption from the mea n intergalac¬ 
tic value, introduced bv IProchaska etld] (l2013ll . This quan¬ 
tity references the equivalent width to the baseline value 
expected from the IGM, 


Spibx, Av) 


W{b±, Av) - WIGM 
AA — WIGM 


(3) 
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Here AA = XaAv/c and wigm = AA[1 — exp(—reff)], where 
Teff is the e ffective optica l dept h of the IGM. We adopt the 
values from I Becker et aP (l2013l i for Teff. Our simulations re¬ 
cover the mean transmission according to these valnes to 
better than 2-4 per cent at 2.0 < 2 < 2.4, and to better 
than 1 per cent at 2.4 < z < 2.8. Rather than renormaliz¬ 
ing the photo-ionization rates in each simulation to recover 
the values of t^s exactly, we take these small deviations to 
be the uncertainty in the model predictions of the amount 
of IGM absorption. All the error bars shown for the model 
predictions are errors in the mean. The error in the mean 
IGM effective optical depth may be added linearly to the 
error bars as a systematic nncertainty. 

Lastly, in the Appendix we present a series of conver¬ 
gence and parameter tests on these statistical measures. In 
summary, we find: (a) a box size of 60 Mpc is preferred over 
30 Mpc for ENZO simulations to capture statistics on haloes 
with Mf, > 10 ^^ Mq; (b) at the resolution of the simulations, 
radiative transfer effects negligibly affect the predicted inte¬ 
grated Hi absorption statistics; (c) when excluding winds, 
the GADGET-3 results for 5f are little changed if the quick 
Lya method is used to convert gas to st ars rather than the 
algorithm of lSpringel fc HernauistI (l2003l L (d) for uncertain¬ 
ties au in galaxy systemic velocities and a velocity window 
An, 5f is insensitive to an for values gh < Au/2; and (e) 
for Oh ~ 500 km s“^ , a velocity window An = 2000 km s“^ 
provides a good compromise between a sensitive dependence 
of on impact parameter and small variance in its values. 


3 RESULTS 

3.1 An illustrative example of Lyo absorption 
around a halo 

It is instructive to first consider the Lyo absorption around 
an individual halo in the cosmological simulations. The gas 
velocity field near a galactic halo is in general complex, in¬ 
volving Hubble flow on large scales, and on smaller, cos¬ 
mological infall, shocks and streams of merging material. It 
is important to include a sufhciently large volume around 
a halo to capture these effects on the simulated spectra. 
Since the spectra are measured in velocity space, large pe¬ 
culiar motions will furthermore scramble spatial and veloc¬ 
ity information. Particularly large displacements in velocity 
space may result from the cosmological infall around massive 
haloes. 

We illustrate these effects for a 3 x 10^^ Mq halo at 
z — 2 from the 30 Mpc GADGET-3 simulation without a 
wind (GSOsfnw), shown in Fig. [T] Three representative lines 
of sight are shown in the lower panels, one through the 
halo centre-of-mass and two offset transversely by ±1.5 
comoving Mpc. Despite the high halo gas temperatures, 
reaching 10^ K, all three spectra show absorption within 
An ~ ±500 km s“^ of the halo centre-of-mass. This absorp¬ 
tion arises in part as a result of the complex velocity field 
of the gas. The panel second from the bottom in particular 
shows the wide range of positions that give rise to line-of- 
sight velocities matching the systemic velocity of the halo. 
These lines of sight include the region within the turnaround 
radius of the gas, r < rt.a. — 6r„, which is sensitive to the 
star formation and feedback implementation (Paper I); the 
virial radius of the halo in Fig.[T]is r„ = 0.48 comoving Mpc. 


Examining individual absorption features, it is appar¬ 
ent that gas extending out to the turnaround radius con¬ 
tributes to the broad feature at —50 km s“^ < v — Uh.ios < 
150 km s“^ along the line of sight passing directly through 
the halo centre. Gas within the turnaround radius also pro¬ 
duces the weaker feature at v — Uh.ios — 25 km s“^ along 
the line of sight displaced laterally by -1.5 Mpc. The line 
of sight displaced by ±1.5 Mpc instead probes part of the 
cosmic web. The very broad absorption feature covering 
—550 km s“^ < V — Uh.ios < —250 km s“^ originates in the 
broad filament attached to the top left of the halo, extend¬ 
ing over —8.5 Mpc < Rios < —3.9 Mpc, as shown in the 
panel second from the bottom. It is dominated by a com¬ 
plex region 30 kpc across (proper) with a velocity width 
of about 160 km s“^ . The resulting absorption feature has 
an equivalent width of 1.2 A and an Hi colnmn density of 
Ahi = 8 X 10^^cm“^, form ing a Lyman limit system along 


this line of sight (see also iFaucher-Giguere fc Kere^ 12011 
Fumagalli et al.ll201ll : | van de Voort et al. I I2OI2I : IShen et all 


2013lf . fAlloi^g for radiative transfer effects, as discussed ir 


the Appendix, increases Nm to 4 x 10^® cm but has little 
effect on the spectrum because the feature is saturated.) 


3.2 Median Lya optical depth 

We now turn to consider the observational signature of the 
neutral hydrogen around galaxy haloes. The most direct 
measurement of Lyo absorption from the gas around galax¬ 
ies is the optical depth To,. Because of its large variance, the 
median optical depth is a more stable statistic than the av¬ 
erage. Fig.[2]shows the median optical depths, as a function 
of separation from the halo centres, averaged over simulated 
spectra in the redshift range 2.0 < 2 < 2.6, corresponding 
to most of the for eground galaxy redshifts in the sample of 
iRakic et al.l (l2012l l. The left panel averages over line-of-sight 
Hnbble flow distances, defined as rios = vios/H{z), within 
±0.4 comoving Mpc of the halo centres. The right panel av¬ 
erages over impact parameters within ±0.5 comoving Mpc 
of the halo centres. Both panels show a general rise toward 
the halo centres, although with differences in the trends. 
In assigning halo velocities, a random Gaussian error with 
standard deviation of au — 130 km s“^ is included to match 
the data ilRakic et al.ll201^ . 

For comoving distances outside 0.7 Mpc, all the model 
predictions in the left panel are in good agreement, with lit¬ 
tle difference with halo mass. At smaller displacements, the 
GADGET-3 simulation with star formation (G30sfnw) shows 
a systematically larger spread in the median optical depths 
compared with the ENZO simulation (E60_1024). The absorp¬ 
tion increases approximately monotonically with halo mass. 
Genera lly good agreement is found in comparison with the 
data of iR.akic et al.l (|2012l l. except for the innermost point, 
where the observations show higher absorption than the sim¬ 
ulations. 

The right panel shows a plateau in the absorption for 
Hubble flow distances within ±1 Mpc of the halo centres, 
and a decline beyond. As in the left panel, the median optical 
depths in the GADGET-3 simulations exceed those of the ENZO 
simulation, with the 10^^ ® M© haloes with winds producing 
the greatest absorption. Outside of 3 Mpc, the median opti¬ 
cal depths are in good agreement b etween all the mod els, as 
well as with the measurements of iRakic et al.l ll2012l L The 
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Figure 1. Top panel: Hi-weighted temperature and proper velocity fields around a 3 x 10^^ Mq halo at 2 = 2 from the GADGET-3 30 Mpc 
box simulation without winds (G30sfnw). The velocities are shown relative to the centre-of-mass of the halo, where an arrow length 
equivalent to the spacing of the axis tic marks corresponds to 1000 km s“^ . The halo virial radius is represented by a circle. The three 
Lyo; absorption spectra shown in the lower panels are drawn along the lines running parallel to the x-axis. Lower panels: line-of-sight 
plots showing the Hi fraction, Phi/{Ph}» gas temperature T, as functions of the line-of-sight velocity offset Au^og from the halo 
centre-of-mass at jos- In the second panel from the bottom, Auios is shown as a function of the comoving distance, R\os^ from the halo 
centre-of-mass. The shaded region indicates the size of the virialized zone. Note that regions at multiple values of -Rios have Auiqs — 0? 
such that gas outside the hot central region contributes to absorption near and across the systemic velocity of the halo. The bottom panel 
shows the corresponding spectra. The open circles and solid black lines are for a line of sight passing through the halo centre-of-mass. 
The inverted cyan triangles (upright blue triangles) and dashed cyan (dot-dashed blue) lines are for a line of sight laterally offset from 
the halo centre-of-mass by -1-1.5 ( —1.5) comoving Mpc. 
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Figure 2. Left panel: median Lya optical depth, averaged over line-of-sight velocities relative to the halo systemic velocity corresponding 
to Hubble distances of up to 0.4 Mpc (comoving), against projected separation from the halo. The results are averaged over all haloes 
in a 0.2 dex wide mass bin centred at the indicated mass. The model predictions agree with the data except within the typical virial 
radius, ry ~ 150-330 kpc (comoving). Right panel: median Lya optical depth, averaged over projected separations up to 0.5 Mpc 
(comoving), against line-of-sight Hubble displacement from the halo. Except for the regions along the spectra corresponding to Hubble 
distances beyond ±3 Mpc (comoving), the models underpredict the measured absorption. In both panels, halo velocities include a random 
component drawn from a Gaussian distribution with standard deviation cth = 130 km s“^ . The data points (large filled circles; red) are 
from Rakic et al. (2012). Note that for clarity of presentation, here and for other figures in this paper, the points in each mass bin are 
grouped with slight offsets using cubic spline interpolation. 
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measured values further in, however, greatly exceed all the 
model predictions. 

Similar agreement with the data for the m edian optical 
as a fu nction of impact parameter was found bv iRakic et al.l 
ll2013l) . with their hydrodynamical simulations also falling 
well short of the measured median optical depths as a func¬ 
tion of Hubble separation along the lines of sight. They point 
out that the discrepancy is difficult to assess statistically be¬ 
cause the data are highly correlated along the line of sight. 
Agreement improves when averaging over increasingly larger 
impact parameters. Rather than pursuing further compar¬ 
isons with median optical depth measurements, however, we 
turn next to integrated spectral quantities across wider line- 
of-sight velocity intervals. We find that these provide a more 
direct means of discriminating between models. 

3.3 The velocity integrated Lya equivalent width 

I Rakic et al.l (l2012l ) provide absorption equivalent widths 
within 1000 km s“^ wide velocity windows centred on the 
systemic velocities of the galaxies in their sample. In Fig. [3l 
the equivalent widths (see Eq. O within a spectral window of 
width Av = 1000 km s“^ centred on the halo velocities from 
the 60 Mpc ENZO simulation (E60_1024) and the GADGET-3 
star formation simulations both without (G30sfnw) and with 
(GSOsfw) supernova feedback are shown as a function of im¬ 
pact parameter b±, averaged over all haloes in halo mass 
bins centred at logjQ Mji = 11.1, 11.5 and 11.9 of width 


A logjQ M — 0.2. The simulated spectra have been ad¬ 
justed by applying a correction factor of 0.804 to the fluxes . 
This matches the correction applied bv iRakic et al.l (|2012| ) 
to allow for errors in t he continuum level in the earlier, 
lower resolution data of ISteidel et al.l (l2010l ). The equiva¬ 
lent width values rise towards smaller impact parameters. 
At b± > 0.4 Mpc (comoving) there is little difference be¬ 
tween the model predictions. At all impact parameters, the 
ENZO predictions are nearly independent of halo mass, with 
the results for the 10^^'^ Mq mass bin slightly smaller than 
the others. By contrast, the GADGET-3 models show a sys¬ 
tematic increase with halo mass for Mu > 10^^’^ M© within 
the inner 0.4 Mpc. For reference, the virial radius ranges over 
150 — 330 kpc for 10^^ — 10^^ Mq haloes0The models includ¬ 
ing supernova feedback predict somewhat larger equivalent 
widths compared with the same mass halo with star forma¬ 
tion alone. This produces a degeneracy in the predictions: a 
given equivalent width value may be produced either by a 
halo with star formation and no wind or a lower mass halo 
with star formation and a wind. 

Gomparison wi th the equivalent width mea surements of 
ISteidel et al.l ll2010l) (large red filled circles) and iRakic et al.l 


^ The comoving virial radius is Vy ~ 0.33M^^^^2 Mpc for a halo 
mass 12 Mq. The turnaround radius, where the infall pe¬ 

culiar velocity cancels the Hubble flow, is at rt a — ^'f'v (Paper 

I)- 
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bj^(conioving Mpc) 


Figure 3. Rest-frame equivalent width (A) within a velocity 
window Av = 1000 km centred on the halo centre-of-mass 
velocity as function of line-of-sight impact parameter fey for halo 
masses logj^Q Mj, = 11.1 (triangles; blue), 11.5 (squares; magenta) 
and 11.9 (circles; black) at redshifts 2.2 < z < 2.5. The data are 
from Steidel et al. (2010) (large solid circles; red) and Rakic et 
al. (2012) (large filled squares; green). 


(|2012|) (large green filled squares) shows excellent agreement 
with the model predictions outside the virial radius. The 
differences between the feedback and non-feedback models 
are too small for the da ta to discriminate between them. 
iTrainor fc SteiddI (l2012l) estimate a median halo mass of 
10^^'® Mq for the survey, consistent with the findings here. 

3.4 Deviation of Lya absorption from the mean 
IGM 

The absorption excess Sp relative to the mean IGM, as quan¬ 
tified by the effective optical depth Teff, averaged over a 
spectral window Av = 1000 km s“^ at 2 = 2.2, is shown 
in Fig. |4] Results are displayed for the ENZO 60 Mpc box 
(E60_1024) and the GADGET-3 simulation with star forma¬ 
tion but no winds (GSOsfnw), averaged over all haloes in 
mass bins centred at logj^g Mh = 11.1, 11.5 and 11.9 of width 
A logj^g M = 0.2. At projected separations b± > 300 kpc 
(comoving), all the models agree, showing a rise towards 
the halo centres (left panel). The non-zero values rising 
towards the halo centres demonstrate excess absorption 
above the average IGM value in an extended region outside 
the haloes, well outside the turnaround radii of the haloes 
(rt.a. ~ 1.8 Mpc comoving for Mh = 10^^ ® Mq). While the 
results for the ENZO simulation show essentially no variation 
with halo mass, the absorption systematically rises within 
the virial radius with halo mass in the GADGET-3 simula¬ 
tions for haloes more massive than 10^^'^ Mq. At projected 
separations smaller than the virial radius, the absorption in 
the GADGET-3 simulation lies below the ENZO values for halo 


masses below 10^^'® Mq. We note that star formation in the 
GADGET-3 simulation substantially reduces the gas density 
within the virial radius compared with the gas-conserving 
ENZO simulations (Paper I). 

The excesses increase slowly with redshift, as shown in 
the right panel of Fig. [J] rising by only about 50 per cent 
from 2 = 2 to 3. The trend of increasing absorption with halo 
mass near the halo centres weakens for the GADGET-3 models 
with increasing redshift, except that the absorption for the 
10^^ ® M© mass bin tends to stay the largest. The different 
models also continue to show the same level of absorption 
at large distances, here shown by the 10^^ ® M© haloes at 
b± — 1.5Mpc. 

Thus far we have compared the simulations with Lya 
absorption data around star-forming galaxies. In the left 
panel of Fig. [5] we instead compare wi th the 5f mea¬ 
surem ents surrounding QSOs reported by IProchaska et al.l 
ll2013t) at 2 Rs 2.4, with a velocity window Av = 
2000 km s~^ . A random component drawn from a Gaussian 
distribution with an = 520 km s“^ is added to the halo ve¬ 
locities to match the typical errors in the measured halo 
redshifts. The predictions for the Mh > 10*^® M© haloes in 
the ENZO simulation lie lower than those for the GADGET-3 
simulations within the turnaround radius, rt.a. ~ 1-8 Mpc, 
although the smaller box for the GADGET-3 simulations re¬ 
sults in larger uncertainties. The measurements outside the 
virial radius agree with all the models. The innermost point, 
however, is consistent only with the Mh > 10^®M© halo pre¬ 
dictions from the GADGET-3 simulations, although the uncer¬ 
tainties are large because of the small number of haloes. 

The fluctuations in the mean transmitted flux, given by 
(Tp = {{F — (F))®)^/®, are shown in the right panel of Fig. [S] 
at 2 = 2.4. The ratio ap /(F) is nearly constant at 0.10 — 0.15 
with bp for all the models. This is substantially smaller than 
the measured fluctuations, especially within the turnaround 
radius. Although the reported values include continuum er¬ 
rors, so that they may perhaps be conservatively taken as 
upper limits, the discrepancy may indicate the need for ad¬ 
ditional astrophysical effects which increase the variance in 
transmitted flux between lines of sight, such as QSO beam¬ 
ing. 

Finally, a further independent test of the models is 
provided by estimating the covering fractions fc for high 
equivalent width absorption. We define an excess integrated 
equivalent width relative to the mean IGM as 

\ /■t'halo-l-AD/2 

Awo{bp, Av) = — / dv e — e ^ (^4) 

where t^h is the IGM effective optical depth. The covering 
fraction is then the ratio of the number of lines of sight with 
Amo above a given threshold to the total number of lines of 
sight. 

The results for Amo > 1.0 A are shown in Fig. |6]for the 
ENZO 60 Mpc box (E60_1024) and the GADGET-3 simulations 
with (G30sfw) and without (G30sfnw) winds. The left panel 
is for absorption in the velocity window Av = 1000 km s“^ 
at 2.0 < 2 < 2.6, and the right for Av = 2000 km s“^ at 
2 = 2.4. The covering fractions rise toward the halo centres 
in all the models. Within the virial radius, the GADGET-3 
models show an increase with halo mass and generally sys¬ 
tematically higher values for the models including wind feed¬ 
back, especially for the narrower velocity window in the left 
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Figure 4. Fractional absorption excess Sp relative to the mean IGM absorption, for a spectral window Av = 1000 km s~^ wide across 
the halo systemic redshifts for the ENZO 60 Mpc box simulation (E60_1024) and the GADGET-3 simulation with star formation but no 
wind (GSOsfnw). Left panel: The variation of 5p with projected impact parameter bj_ and halo mass at z = 2.2. Absorption exceeds 
the average IGM contribution out to several comoving Mpc. Right: The evolution of Sp with redshift for selected values of b± and halo 
mass. The symbols are as indicated on the diagram. In both panels, halo velocities include a random component drawn from a Gaussian 
distribution with standard deviation cjh = 130km s~^ . 
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Figure 5. Left panel: Fractional absorption excess Sp relative to the mean IGM absorption, for a spectral window Av = 2000 km s“^ 
across the halo systemic velocities for the ENZO 60 Mpc box simulation (E60_1024) and the GADGET-3 si mulations with (G30sfw ) and 
without (G30sfnw) a wind. The data points in the left panel are for absorption around QSOs taken from [^ochaska et al.l ll2013lj . with 
the error bars showing the errors in the mean. Right panel: The fluctuations relative to the mean transmitted flux, ap / (F) (see text 
for details). Both Sp and ap/{F) are shown against projected impact parameter b± at z = 2.4. The halo velocities include a random 
component drawn from a Gaussian distribution with standard deviation cth = 520 km s~^ . 
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with standard deviation cth = 130 km s~^ (left panel) or cth = 520 km s“^ (right panel). 
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Figure 7. Covering fraction fc for discrete absorption systems with column densities A^hi > 10^® cm“^ within a spectral window of 
width Av = 600 km s“^ across the halo systemic velocities (left panel) and A^hi > 10^^’^ cm~^ for a window Av = 3000 km s“^ (right 
panel) for the GADGET-3 simulations with (G30sfw) and without (G30sfnw) a wind. The covering fraction is shown against projected 
impact parameter bj^ at z = 2.3. The halo velocities include a random component drawn from a Gaussian distribution with standard 
deviation cth = 130 km s“^ (left panel) or cth = 520 km s“^ (right panel). 
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panel. The more massive haloes with winds achieve covering 
fractions fc > 0.5 well within their virial radii. 

While the integrated excess equivalent widths are 
straightforward to measure from observations, the focus of 
recent literature has been the covering fractions of identi¬ 
fied discrete absorption systems. A direct comparison with 
the simulations would require absorption line identifica¬ 
tion and fitting, which we defer to later work. Here, we 
use the high spatial resolution of the GADGET-3 simulations 
to provide illustrative comparisons by identifying discrete 
contiguous Hi systems along the lines of sight. In Fig. 0 
we compare the predicted cover ing fractions with the ob¬ 
servations of iRudie et al.l ll2012ll for systems with H i col¬ 
umn densities Ahi > 10^® cm“^ within a velocity window 
An = 600 km s“^ wide centred on the gal axy systemic ve¬ 
locitie s (left panel), and with the data of IProchaska et aP 
ll2013ll for “optically thick” absorbers in a velocity win¬ 
dow An = 3000 km s“^ wide centred on the galaxies (right 
panel). These authors define optically thick absorbers vari¬ 
ously as those showing obvious damping wings, Lyman limit 
absorption, strong low-ionization metal absorption or (if not 
classifiable otherwise) exhibiting a single strong Lya feature 
with Wo > 1.8 A. We select discrete systems from the simu¬ 
lations with column densities Ahi > 10^^'® cm“^. 

The covering fraction of systems with Ahi > 10^® cm~^ 
is small compared with the measurements for most of the 
models, particularly within the virial radius. The excep¬ 
tion is for massive haloes. Mu > 10^^'® M©, for which the 
model with a wind shows a rapid rise in the covering frac¬ 
tion toward the halo centre close to the measured values. 
The covering fractions for absorbers optically thick at the 
Lyman limit (Ahi > 10^^'® cm~^) l i e syst ematically below 
the observations of IProchaska et ^ (l2013l L While the cov¬ 
ering fraction for Mh > 10^^ Mq approaches the measured 
values within the spread from the simulation, the uncer¬ 
tainties in the mean are large due to the small number 
of haloes available in this mass bin (a total of 12). This, 
along with the heterogeneous definition of the observed op¬ 
tically thick absorbers, makes it difficult t o quantify statis¬ 
tically the level of disagreement. Moreover. IProchaska et al] 
(l2013ti caution against errors due to continuum placement 
and line-blending. Similarly low c overing fractions com ¬ 
pare d with observations are f ound bv iFumagalli et al] (l2014l f 
and iFaucher-Giguere et al.l (l2014h for halo masses exceed¬ 
ing 10 ^^ Mq. By contrast, a recent study bv iRahmati et aP 
(I2OI5I I finds close agreement with the measured values, 
which they attribute to their feedback model. Their compar¬ 
ison, however, is for larger mass haloes, Mh > 10^^'® Mq, 
captured within their larger simulation volume. It should 
also be noted that the simulations other than ours compute 
the covering fractions based on integrated line-of-sight col¬ 
umn densities rather than discrete absorption systems. 


4 DISCUSSION 

4.1 Radial characteristics of the absorption 

We have compared the statistics for Lya absorption along 
lines of sight passing through the environments of galaxy 
haloes using three simulations: an ENZO simulation without 
star formation and two GADGET-3 simulations with star for¬ 
mation, one of which also allows for feedback in the form of 


supernovae-driven winds. We found in Paper I that the gas 
properties of the haloes largely agree beyond the turnaround 
radius, regardless of the inclusion of star formation or winds. 
The results presented here confirm that the predicted ab¬ 
sorption line statistics for lines of sight passing outside the 
turnaround radius are similar. 

Co mparisons w i th th e Lyg absorption statistics in 
SFGs dS teidel et al.l 201 ol : iRakic et al.l I 2 OI 2 II and QSOs 
(IProchaska et al.l l2013ll . demonstrate the simulations suc¬ 
cessfully recover the absorption statistics outside the halo 
virialized regions. Within the virial radius, however, the 
measured amount of integrated absorption rises to a level 
none of the models reproduce. A similar trend is found us¬ 
ing the median Lya optical depth data for SFGs (Fig. [2j. 

The various absorption line statistics nevertheless show 
some common trends among the models. The predictions of 
all the models, whether or not they include star formation or 
wind feedback, all agree, within the error bars, outside the 
turnaround radii of the haloes. The models also agree well 
with the data for both SFGs and QSOs in this region. As 
the virial radius is approached, however, the models begin 
to diverge. 

The gas around galaxies has generally been divided 
between a circumgalactic medium wit hin a distance of ~ 
300 kpc (proper) fro m a galaxy (e.g. ISteidel et al.l I 2 OIC 1 I : 
IProchaska et al.ll2013l L and the IGM on larger scales. Based 
on our simulation results, we rather suggest it is useful to 
consider three distinct regions: the inner virialized region, 
an intermediate region we call the ‘mesogalactic medium’ 
(MGM), extending between the virial radius and the IGM 
(r„ < r < 12r„), and the IGM well outside galactic haloes 
(r > 12r„). The virialized region is distinguished as the 
active theatre within which star formation and feedback 
most affect the hydrogen absorption signatures. Here, dis¬ 
crete optically thick absorption systems, especially damped 
Lya absorbers, may make substantial contributions to the 
integrated amount of absorption. The MGM is the region 
around galaxies exhibiting enhanced absorption over the dif¬ 
fuse IGM, but little affected by star formation or feedback. 
It originates from gas in the extended haloes around galax¬ 
ies (Paper I), before merging into the diffuse IGM on larger 
scales. 


4.2 From where does the measured absorption 
arise? 

Many of the Hi absorption observational studies around 
galaxies have focussed on the search for evidence of the cold 
streams predicted to penetrate into the virialized regions of 
moderate mass galactic haloes [Mh < 10^^ Mq). Establish¬ 
ing that the absorption arises from inward streaming gas 
is observationally challe nging; indeed, the measured kine¬ 
matics favour outflows (ISteidel et al.l 12010 1. The measured 
covering fractions of cold {T « K) ga s appear to exceed 
theoretical expectations in both S FGs (iRakic et al.l I2OI2II 
and QSOs (IProchaska et al]|2013ll . The excess absorption 
in QSOs is especially intriguing, as the gas interior to the 
virial radius is expected to be shock heated to such high 
temperatures as to be in collisional ionization equilibrium, 
with greatly reduced levels of neutral hydrogen. As shown in 
Fig-Ill fhe absorption signal as a function of impact param¬ 
eter is complicated by the complex peculiar velocity field of 
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the gas. The large peculiar motions of the gas in the vicin¬ 
ity of galaxies can shift the absorption contribution of gas 
parcels distant from the galaxies into and out of the velocity 
window used to measure the absorption signature. 

In order to gain some insight into the origin of the mea¬ 
sured absorption signatures, we construct spectra from the 
ENZO 60 Mpc simulation after masking all the Hi within the 
virial radius of haloes more massive than Mh > 10^^ Mq. 
The resulting effect on the equivalent width in a velocity 
window An = 1000 km s“^ wide centred on the galaxy ve¬ 
locities is shown in the left panel of Fig. [8] (open symbols). 
While the contribution of gas external to the virialized re¬ 
gion to the absorption signature decreases towards small im¬ 
pact parameters for b± < r„, it never vanishes. At least 70 
per cent of the absorption at << r„ arises from gas out¬ 
side the virialized zone. 

In terms of the absolute equivalent width values, Fig. [3] 
shows that within the virial radius the simulations under¬ 
predict the measured equivalent width by Awo — 1 A com¬ 
pared with observations. About half the mean absorption 
is thus missing within the virial radius. Zoom-in simula¬ 
tions, achieving a baryon mass resolution of about 10® Mq 
for haloes more massive than 10^^ Mq, find covering frac¬ 
tions for_saturated_^sorgtionjlnessmallemth^ 50 pe r cent 
(e.g. iFumaealli et al.|[201ll . l2014l : iGoerdt et al.l 12012 1. Un¬ 
til damping wings appear (Nm > 10^^cm“^), an individ¬ 
ual saturated line will contribute only about Awo — 0.5 A. 
Given the covering fraction, this is too small to make 
up the difference. Systems with column densities Nm > 
10^® cm“^, or equivalent widths of 2 — SA, are rare, with 
covering fractions well under 10 per cent in these simu¬ 
lations, so that including systems showing damping wings 
c annot fully make up the differ ence either. The simulations 
of iFaucher-Giguere et al.l ll20l4 l. with a baryon resolution as 
small as 4 X10^ Mq , however, show a higher incidence rate of 
damped absorbers. For haloes more massive than 10^^ Mq, 
the covering fractions for saturated lines reach 60-80 per 
cent, and 20-40 per cent for systems with A^hi > 10^® cm~^. 
This may just be sufficient to make up the missing absorp¬ 
tion from the virialized zone. If the relatively rare systems 
with damped wings account for the missing absorption, then 
wide variations in the integrated absorption between differ¬ 
ent lines of sight is expected, with the absorption found in 
the simulations presented here providing a more constant 
baseline level. 

Masking out the region within the turnaround radius 
(r < 6r„) further reduces the signal, but the contribution 
to absorption within b± < 6r„ from material outside the 
turnaround radius is still non-negligible. Gas at r > 6r„ 
contributes at least 25 per cent of the absorption signal we 
obtain even at << r„ (by comparison, less than 10 per 
cent of the signal derives from the diffuse IGM, at r > 12r„). 
At transverse separations outside the turnaround radius, the 
full absorption signal is still not achieved, demonstrating 
that about 20 per cent of the absorption along a line of 
sight passing through the mesogalactic region of a given halo 
arises from mesogalactic gas in other haloes. For a halo of 
mass Mh = 10^^'® Mq, this corresponds to absorption from 
material at a distance beyond 1.8 Mpc (comoving), or a ve¬ 
locity offset of 130 km s“^ , well within the velocity window 
Av = 1000 km s“^. 

Increasing the velocity window to Av = 2000 km s“^ 


(Fig. [8] right panel) produces similar trends. Results for 
halo masses corresponding to the host galaxies of QSOs 
(12.0 < logjQ Mh < 12.6) are now also shown. For the in¬ 
creased velocity window, even for these massive haloes 80 
per cent of the absorption signal at small transverse impact 
parameters {b± « r„), originates in gas outside the virial¬ 
ized region. About half the signal for lines of sight passing 
within the turnaround radius originates from gas beyond the 
turnaround radius, while 20 per cent of the absorption for 
lines of sight passing through the virialized and mesogalactic 
regions originates from gas in the diffuse IGM. 


5 CONCLUSIONS 


We perform large-scale cosmological hydrodynamical 
simulations using two numerical schemes, GADGET-3, 
including star formation using the prescription of 
ISpringel fc Hernauisd ll2003l l. with and without supernovae- 
driven wind feedback, and ENZO without star formation 
as a control case. The simulations allow us to assess the 
impact of star formation and wind feedback separately on 
the Hi absorption statistics around galaxies and QSOs. The 
statistical measures we consider are based primarily on the 
integrated absorption properties within velocity windows 
centred on the systemic velocities of the galactic haloes: 
equivalent width values, fractional absorption excesses 
relative to the mean IGM and fluctuations in the absorbed 
flux. We also provide statistical predictions for covering 
fractions of integrated equivalent width values over the 
velocity windows and for some discrete absorption systems. 
We compare th e simu l ation results with t he m easurements 
ofISteidel et al.l (|201(t l. iRakic et aP ll2ni2h andlRud e et alJ 
(|2012ll for star-forming galaxies and IProchaska et al. ( 2013II 
for QSOs over the redshift range 2 < 2 < 3. 

On the basis of the simulated absorption properties, 
we identify three regions in the gas surrounding a galaxy 
with distinct absorption properties, the inner virialized re¬ 
gion, a mesogalactic zone extending from the virial radius 
to twice the turnaround radius (rt.a. — 6r„ « 1 —3 cMpc for 
IqH.o _ ]y[g haloes), and the IGM beyond. The simu¬ 
lations show star formation and wind feedback play only a 
secondary role in establishing the integrated Hi absorption 
signatures compared with the general trend of increasing 
absorption for decreasing impact parameter. All the simu¬ 
lations reproduce the absorption measurements within the 
mesogalactic medium and beyond. This is a non-trivial re¬ 
sult since it extends the success of IGM simulations in the 
context of a AGDM cosmology, from the diffuse IGM to the 
extended haloes of galaxies down to their virial radii, with¬ 
out being subject to the uncertainties of subgrid physics. 
As such, the agreement is a confirmation of the AGDM cos¬ 
mological model for structure formation on comoving scales 
down to 0.4 Mpc around galaxies over the redshift range 
2 < 2 < 3. 

The GADGET-3 simulations, both with and without a 
wind, and the ENZO 60 Mpc box simulation produce compa¬ 
rable equivalent width values in Av = 1000 km s“^ windows 
centred on the halo systemic velocities over the full range of 
impact parameters. Allowing for a wind, all the model pre¬ 
dictions for haloes wi th masses Mh > M q agree 


the measurements of ISteidel et al. ll2010ll and Rakic et al.l 
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Figure 8. Fraction of equivalent width arising from gas outside the virialized region, r > r„ (open symbols), outside the turnaround 
radius, r > 6r„ (solid symbols), and outside the mesogalactic region, r > 12r„ (starred symbols), shown against projected comoving 
impact parameter by at 2 = 2.4. The equivalent widths are computed within spectral windows of width Av = 1000 km s“^ (left panel) 
or 2000km s“^ (right panel) across the halo systemic velocities for the ENZO 60 Mpc box (E60_1024). In both panels, most of the 
absorption for lines of sight passing within the virial radius arises from gas outside the virialized region, and about a third from gas 
outside the turnaround radius. The halo velocities include a random component drawn from a Gaussian distribution with standard 
deviation (th = 130 km s“^ (left panel) or 520km s“^ (right panel). 


ll2012l) outside the virial radius. Within the virialized region, 
the GADGET-3 models show a secondary trend of increasing 
absorption with halo mass and enhanced absorption for a 
given halo mass when wind feedback is included, produc¬ 
ing a degeneracy between halo mass and wind feedback, at 
least for the feedback model we adopt. All the models, how¬ 
ever, underpredict the amount of absorption compared with 
observations for lines of sight passing through the virial¬ 
ized regions. High resolution zoom-in simulations including 
radiative transfer suggest Lyman limit systems, especially 
damped Lya absorbers, may account for the remaining ab¬ 
sorption. 

We also compare our models with the median opti¬ 
cal depth measurements of iRakic et al.l ll2012l) . The models 
agree well with the measurements outside the virial radius, 
where all the models show similar values, but not within, 
similar to the integrated equivalent width comparisons. In 
addition to possibly being a consequence of unresolved ab¬ 
sorption systems, the shortfall may also be in part statisti¬ 
cal, a consequence of line-of-sight correlations in the median 
optical depth val ues on small scales (we refer the reader to 
iRakic et al.ll2013l . for a discussion). 

We also compare the model predicti ons for the frac¬ 
tional absorption excess 5p with the data of IProchaska et ^ 
ll2013l) . who report measured values in velocity windows 
Av = 2000 km s“^ centred on QSOs. In the mesogalactic 
region, both the ENZO and GADGET-3 simulation predictions 
are largely insensitive to the halo mass. Agreement with 
the measured values is achieved even for halo masses below 


Mh > ° Mq, the expected lower halo mass of QSO host 

galaxies. The ENZO simulation underpredicts the amount of 
absorption within the virial radius, as found for SFGs. The 
comparison with the GADGET-3 simulations is inconclusive 
because of the low number of massive haloes in the simula¬ 
tion volume. 

We compute covering fractions for integrated equiva¬ 
lent widths within the velocity windows, as these are read¬ 
ily measured from observations without the requirement of 
identifying and fitting discrete absorption systems in the 
spectra. The sensitivity we find to feedback suggests the 
covering fractions may provide useful constraints on feed¬ 
back models. Since much of the literature focusses instead 
on the covering fraction of discrete absorption systems, we 
also consider illustrative comparisons with observations for 
saturated absorption lines. We find good agreement with ob¬ 
servations of SFGs, but underpredict the covering fraction of 
optically thick absorbers in QSOs, with only marginal agree¬ 
ment with the observations. Because of the low number of 
massive haloes in our simulations, however, the statistics 
are too poor to be conclusive. Agreement may require larger 
mass haloes than may be examined with the limited box 
size analysed here. An accurate determination of the contri¬ 
butions from individual Lyman limit systems and damped 
Lya absorbers to the absorption properties within the viri¬ 
alized region may also require the increased resolution of 
zoom-in simulations. Recovering the full amount of absorp¬ 
tion from virialized gas will likely require additional phys¬ 
ical effects, including alternative sub-grid feedback models. 
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self-consistent radiative hydrodynamics to account for the 
physical response of systems optically thick to ionizing radi¬ 
ation, or p ossibly pressure resulting from a fluctuating mag- 
netic field (IPandev fc Sethill2013l : Chongchitnan fc Meiksiiil 
|2014| . 
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Figure Al. Fractional absorption excess Sp relative to the mean 
IGM absorption for a velocity window An = 2000 km s~^ around 
the halo centre-of-mass velocity. The data are shown as a func¬ 
tion of the line-of-sight impact parameter bj_ for halo masses 
11.4 < logjQ M/i = 11.6 (squares) and 12.0 < log^Q < 12.6 

(circles) at 2 : = 2.3, for the ENZO simulations in boxes of side 
30 Mpc (E30_1024, open symbols) and 60 Mpc (E60_1024, filled 
symbols). While the predictions converge well for moderate mass 
haloes, differences arise for the most massive haloes, although 
the uncertainties are large because of the small number of haloes 
in the 30 Mpc box. Halo velocities include a random compo¬ 
nent drawn from a Gaussian distribution with standard deviation 
(Th = 520 km s“^ . 
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Figure A2. Fractional absorption excess Sp relative to the mean 
IGM absorption within a Av = 1000 km s“^ window centred on 
the halo centre-of-mass velocity. The data are shown as a func¬ 
tion of the line-of-sight impact parameter bj_ for halo masses 
\ogiQ Mh = 11.1 (triangles; blue), 11.5 (squares; magenta) and 
11.9 (circles; black) at redshift z = for the GADGET-3 non-wind 
simulations using the default star-formation prescription of SH03 
(G39sfnw, filled symbols) or quick Lya (G30qLya open symbols). 
The choice of gas removal algorithm has little impact on the 
amount of absorption. Halo velocities include a random compo¬ 
nent drawn from a Gaussian distribution with standard deviation 
(Th = 130 km s“^ . 
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APPENDIX A: CONVERGENCE TESTS ON Hi 
STATISTICS 

The convergence of the fractional absorption excess 5p with 
box size is tested in Fig. EH using the ENZO simulations for 
a velocity window Av = 2000 km and a halo redshift 


uncertainty an = 520 k m s ^ , similar to the observations of 
IProchaska et al.l ll2013ll . Only small changes in (5 f, less than 
0.02, are found for the gas around ~ 10^^ ® Mq haloes. In 
the high mass bin there is considerable Poisson scatter in 
Sf from the smaller box, in which there are only 14 haloes 
with Mh > 10^^ Mq. 

We examine the possible role radiative transfer may 
have on the absorption features using the simplifying ap¬ 
proximation of an attenuated radiation field within systems 
sufficiently dense to be self-shielding to photoionizing rad i- 
ation. We adopt the prescription of iRahmati et all ll2013li . 
using a characteristic self-shielding total hydrogen density of 
0.0064r4 '^^ cm~® for temperature T 4 = T/IO"* K. The effect 
on the mean values of Sf are under 5 per cent, and generally 
less than one per cent. We therefore neglect the effects of ra¬ 
diative transfer in this paper on the integrated amount of ab¬ 
sorption. Radiative transfer does affect the column densities 
of discrete absorption systems, but acts principally by cre¬ 
ating systems with Hi column densities sufficiently large to 
produce damping wings (Nm > 10^® cm“®). Since we cannot 
include the physical response of the gas due to the resulting 
temperature and consequent internal pressure changes these 
systems would undergo, we do not consider them except for 
their small contribution to the lower thresholds of systems 
we do examine, Nm > 10^® cm“® and Nm > cm~® 
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Figure A3. Fractional absorption excess 5f relative to the mean 
IGM for halo masses 12.0 < log^^Q < 12.6 within velocity 
windows of width Ai; = 1000km s“^ (triangles), 2000 km s~^ 
(squares) and 3000km s“^ (circles). The results are centred on 
the halo centre-of-mass velocity as a function of the line-of- 
sight impact parameter b± for the GADGET-3 non-wind simula¬ 
tions (GSOsfnw) at 2 = 2.2. While increasing the velocity window 
reduces the variance, it also suppresses the signal. The window 
Av = 2000 km s“^ is a good compromise. Halo velocities include 
a random component drawn from a Gaussian distribution with 
standard deviation cth = 520 km s~^ . 


Figure A4. Relative flux fluctuation ap/^F) averaged over a 
velocity window Au = 2000 km s“^ centred on the halo centre- 
of-mass velocity for halo masses log^^Q My^ = 11.1 (triangles; blue), 
11.5 (squares; magenta) and 11.9 (circles; black) at b± = 3.8 Mpc 
(comoving). The results are displayed as a function of redshift for 
the GADGET-3 non-wind simulation (G30sfnw). The dashed line is 
the fit from Prochaska et al. (2013) to IGM measurements. It is 
most closely approached for the lower mass haloes. Halo velocities 
include a random component drawn from a Gaussian distribution 
with standard deviation cth = 520 km s“^ . 


in Fig. [71 for which the effects of radiative transfer are in¬ 
cluded. We have also compared results for Sf using Doppler 
and Voigt profile functions for run GSOsfnw. The differences 
were negligible. We use the Doppler profile function for the 
results presented here since it requires an order of magnitude 
less analysis time. The Voigt function was used, however, 
when including the effects of radiative transfer. 

All the GADGET-3 simulations for the results in 
this paper use the sta r formation prescription of 
ISoringel &; HernauistI (|2003l) . In Fig. IA2l we compare the re¬ 
sults for Sf using the quick Lya method instead. Very little 
difference is found, suggesting the results are robust to the 
method of gas removal. 

Increasing the velocity window width Av suppresses 
the values of Sf at all impact parameters within the 
turnaround radius of massive haloes, while reducing the 
spread in values, as shown in Fig. IA3I Decreasing Av 
from 2000 km s“^ to 1000 km s“^ slightly increases the sig¬ 
nal while producing a somewhat wider spread. A window 
width of Av — 2000 km s“^ produces a good compromise 
between these competing effects. Velocity window widths 
A v = 1000 km s ~^ an d 200 0 km s~^ to compute Sf are used 
by iRakic et al.l (l2012l ~) and IProchaska et al.l (l2013h , respec¬ 
tively, much wider than the respective typical halo velocity 
errors of an = 130 km s“^ and 520 km s“^ . We have con¬ 
firmed that a halo velocity uncertainty as large as half the 


velocity window affects the values of Sf by less than two per 
cent. 

Finally, the relative fluctuations in the Lya flux in a 
given velocity window in random diffuse IGM regions in¬ 
crease with redshift. In Fig. IA4I we show the convergence in 
our simulations to the IGM data measured in a spe ctral win¬ 
dow of Av = 2000 km s“^ llProchaska et al]|2013l l for large 
impact parameters around the haloes. The diffuse IGM value 
is approached most rapidly in the gaseous surroundings of 
the lower mass haloes, for which the mesogalactic region is 
smaller. 
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